Simulating radiation and kinetic processes in relativistic plasmas

Context. Modelling the emission properties of compact high energy sources such as X-ray binaries, AGN or γ-ray bursts represents a complex problem. Contributions of numerous processes participate non linearly to produce the observed spectra: particle-particle, particle-photon and particle-wave interactions. In the past decades, numerical simulations have been widely used to address the key properties of the high energy plasmas present in these sources. Aims. This article presents a code that has been designed to investigate these questions. It includes most of the relevant processes needed to simulate the emission of high energy sources. Methods. This code solves the time-dependent kinetic equations for homogeneous, isotropic distributions of photons, electrons and positrons. No assumption is made on the shape of these distributions. Have been included so far: synchrotron self-absorbed radiation, Compton scattering, pair production/annihilation, e-e and e-p Coulomb collisions and some prescriptions for additional particle heating and acceleration. Results. We also present comparisons with earlier works and some examples to illustrate the code computational capacities. Previous results are reproduced qualitatively but some differences are often found in the details of the particle distribution. As a first application of the code, we investigate acceleration by second order Fermi-like processes and we find that the energy threshold for acceleration has a crucial influence on the particle distribution and the emitted